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ABSTRACT 

We have compiled photometric data from the Wide-field Infrared Survey Explorer 
All Sky Survey and other archival sources for the more than 2200 objects in the original 
McCook & Sion Catalog of Spectroscopically Identified White Dwarfs. We applied 
color-selection criteria to identify 28 targets whose infrared spectral energy distributions 
depart from the expectation for the white dwarf photosphere alone. Seven of these are 
previously known white dwarfs with circumstellar dust disks, five are known central stars 
of planetary nebulae, and six were excluded for being known binaries or having possible 
contamination of their infrared photometry. We fit white dwarf models to the spectral 
energy distributions of the remaining ten targets, and find seven new candidates with 
infrared excess suggesting the presence of a circumstellar dust disk. We compare the 
model dust disk properties for these new candidates with a comprehensive compilation 
of previously published parameters for known white dwarfs with dust disks. It is possible 
that the current census of white dwarfs with dust disks that produce an excess detectable 
at K-band and shorter wavelengths is close to complete for the entire sample of known 
WDs to the detection limits of existing near-IR all-sky surveys. The white dwarf dust 
disk candidates now being found using longer wavelength infrared data are drawn from 
a previously underrepresented region of parameter space, in which the dust disks are 
overall cooler, narrower in radial extent, and/or contain fewer emitting grains. 
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1. Introduction 

Dust disks are common in a wide variety of astrophysical situations, including the central 
engines of quasars and active galactic nuclei (e.g., Rowan- Robinson 1977; Antonucci 1993), the 
precursors of planetary system formation around protostars (e.g., Natta 2008) and post-formation 
debris around young stars (e.g., Aurnann 1985; Chini et al. 1991), and even the recently discov- 
ered largest ring of Saturn (e.g., Verbiscer et al. 2009). The first dust disk around a white dwarf 
(WD), G29-38, was discovered by virtue of its infrared (IR) excess over the WD photosphere 
(Zuckerman & Becklin 1987), although it took another decade to cast aside all doubts that the 
excess was truly due to dust and not an unresolved brown dwarf companion (Koester et al. 1997; 
Kuchner et al. 1998). Debes & Sigurdsson (2002) and Jura (2003) developed a model for the origin 
of WD dust disks involving tidal disruption of a comet or asteroid perturbed into the WD Roche 
lobe due to the gravitational influence of a remnant planetary system containing at least one mas- 
sive planet (Debes et al. 2012). It wasn’t until 2005 that the second WD with a dust disk, GD 362, 
was discovered (Becklin et al. 2005; Kilic et al. 2005). By the end of 2010, 20 dusty WDs were 
known (Table 5.1 in Farihi 2011), largely owing to sensitive IR observations from the Spitzer Space 
Telescope. 

Predating the discovery of dust around WDs, it was known that a small fraction of WDs show 
absorption lines of metals in their optical and UV photospheric spectra (e.g., Lacombe et al. 1983; 
Shipman &; Greenstein 1983; Zeidler-K.T. et al. 1986). These lines originate from “pollution” of a 
WD’s geometrically thin, high density (but non-degenerate) atmosphere, which in most objects is 
otherwise pure hydrogen or helium. Gravitational settling times in hydrogen-rich (DA) WD atmo- 
spheres are very short (a few days to < 1000 yr), so metals quickly diffuse out of the photosphere 
unless replenished. Thus, the observed metals were thought to be supplied by ongoing accretion 
from the ISM (Sion et al. 1990). This explanation was problematic; notably, the required accre- 
tion rate is high (> 10 8 ' 5 g s -1 compared to ~ 10 7 g s^ 1 expected from ISM accretion) and it is 
difficult to explain the relative elemental abundances of the accreted material, which do not match 
equilibrium ISM values (see Farihi et al. 2010 and the review and discussion in Sections 5.2.4 and 
5.6.6 of Farihi 2011). 

It is a testament to the strength of the asteroid disruption model for WD dust disks that it 
also explains the metal-rich WDs, via accretion from circumstellar dust. Zuckerman et al. (2007) 
showed that the relative abundances of accreted metals in the dusty WD GD 362 closely match those 
of the terrestrial planets. Analyses of two metal polluted WDs (GD 61, NLTT 43806) suggest that 
the accreted dust derived from an asteroid whose origin was in the outer layers of a differentiated 
planet, in which the heaviest elements had sunk to the core, leaving a lithosphere rich in Ca and 
possibly water (Farihi et al. 2011; Zuckerman et al. 2011). Thus, observing WDs with dust disks 
is directly linked to determining how the chemical diversity of planetary systems can influence the 
probability that some planets support life. 

In order to facilitate this, a large sample of WDs with dust disks is desired. Consequently, 
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we have been carrying out the WISE InfraRed Excesses around Degenerates (WIRED) Survey, 
utilizing photometry from the Wide-field Infrared Survey Explorer (WISE), a NASA medium class 
Explorer mission launched on 14 Dec 2009 (Wright et al. 2010). WISE mapped the entire sky 
at 3.4, 4.6, 12, and 22 pm (Wl, W2, W3, and W4 bands, respectively) with 5 a point source 
sensitivities of approximately 0.08, 0.11, 1, and 6 mJy, respectively. Complete sky coverage was 
achieved in mid-July 2010. The WIRED Survey has the goals of characterizing WD stars in the 
WISE bands, confirming objects known to have IR excess from past observations (Spitzer, 2MASS, 
UIvIDSS, etc.), and revealing new examples of WDs with IR excess that can be attributed to 
unresolved companions or circumstellar debris disks. We are utilizing target lists drawn from 
cataloged WD samples (e.g., from the Sloan Digital Sky Survey, McCook & Sion, etc.). To date, 
we have published results from WIRED for the SDSS WD sample that have nearly tripled the 
number of known WDs with circumstellar dust disks and increased the number of WD + brown 
dwarf binaries by almost an order of magnitude. We now present initial results from examining 
the McCook & Sion (1999) (henceforth, MS99) catalog of spectroscopically identified WDs, which 
resulted in the new identification of seven WDs with IR excess indicating the likely presence of a 
circumstellar dust disk. 


2. Targets and Data 

MS99 contains 2249 (optical) spectroscopically identified WDs. The updated and online version 
of MS99, the Villanova University White Dwarf Catalog 1 (henceforth, MSonline) currently contains 
over 14,000 entries. For our purposes, we have used only those targets contained in the original print 
publication of MS99 (most of the new WDs listed in MSonline are objects discovered by SDSS and 
are covered in our previous WIRED paper; see Debes et al. 2011). Since the publication of MS99, a 
number of the WDs were subsequently reclassified as non- WDs (e.g., quasars) or nonexistent (e.g., 
some WDs are listed twice in MS99 with different names, such as WD 2321-549 = WD J2324-546), 
leaving 2202 viable targets. To further narrow the target list, we considered only the 1474 WDs 
from MS99 for which Hoard et al. (2007) found 2MASS near-IR detections. 

As a first pass at identifying IR-bright (or otherwise “interesting”) WDs from this input list, 
we utilized the following selection criteria for each target: (1) 2MASS-J, Wl, and W2 photometry 
exists, (2) the WISE color index (Wl— W2) > +0.3 mag, and (3) the signal-to-noise ratios of the 
Wl and W2 detections are both > 7. 

The IR color-color diagram of the MS99 WDs (see Figure 1) demonstrates that the 2nd criterion 
selects the majority of known dusty WDs, while avoiding the bulk of the “uninteresting” WDs. 
The 3rd criterion excludes color-selected targets whose redness is spurious, resulting from low S/N 
photometry. A broad 10 pm silicate emission feature is a hallmark of circumstellar dust around 


1 See http : //www . astronomy . villanova.edu/WDCatalog/ index.html. 
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WDs (Jura et al. 2009) and falls into the W3 band. In principle, this could offer an additional 
selection criterion for identifying WDs with dust. In practice, however, we found that due to the 
lower sensitivity of WISE in the W3 band compared to W1 and W2 (cf. while 67% and 58% of 
our input sample have a W1 and/or W2 detection, respectively, only 16% have a W3 detection), 
almost 92% of the targets that have a cataloged W3 detection are already selected by our 1st 
criterion. Only two of the targets with a W3 detection that were not selected by our 1st criterion 
have S/N > 7 in W3, and both of them are unusable: WD 1919+145 is contaminated by a nearby 
source (Mullally et al. 2007) and WD 2110+300 is in an unresolved binary with a G-type giant 
star (£ Cygni; Griffin & Keenan 1992). In addition, closer inspection of the WISE images of our 
targets shows that as many as ~50% of the cataloged W3 “detections” (especially at low S/N) are 
probably unreliable (e.g., due to bright, structured background or nearby sources) and should be 
treated as upper limits (e.g., see Section 2.1). 

Incidentally, there are some features of the IR color-color diagram that we will not discuss in 
detail, but are worth noting: 

• There are two principal loci of WDs (plotted as small gray circles) : a large one centered around 
(J— Wl) ~ —0.25 mag, (Wl— W2) « 0 mag, and a smaller one at (J— Wl) « +1.0 mag, 
(Wl— W2) ss +0.2 mag. The former is the locus of “naked” WDs, while the latter is the 
locus of unresolved binaries containing a WD and a low mass main sequence star. We found 
a similar distribution of the majority of points in the r, i, J, Wl, W2 color-color planes for the 
targets selected from the SDSS Data Release 7 preliminary WD catalog (Debes et al. 2011). 

• The known WDs with circumstellar dust disks form a broad sequence extending from the locus 
of naked WDs (with blue colors) to the upper right (red) corner of the color-color diagram. 
The new candidate WDs with dust disks reported here (see below) follow this sequence at 
the red end, and broaden it at the blue end. 

Our criteria initially resulted in the selection of 28 WDs from the MS99 sample. Seven 
of these were excluded because they are known WDs with circumstellar dust disks: WD 0408- 
041 (Kilic et al. 2006), WD 0843+516 (Gansicke et al. 2012), WD 1015+161 (Jura et al. 2007a), 
WD 1116+026 (Jura et al. 2007a), WD 1150-153 (Kilic &: Redfield 2007), WD 1541+650 (Kilic et al. 
2012), and WD 1729+371 (Becklin et al. 2005; Kilic et al. 2005). Five WDs were excluded for 
being known central stars of planetary nebulae (CSPNs): WD 0558-756 (Hcnize &; Fairall 1981; 
Rauch & Werner 1993), WD 0950+139 (Liebert et al. 1989), WD 1821+643 (Starrfield et al. 1985), 
WD 1958+015 (Napiwotzki & Schoenberner 1995), and WD 2333+301 (Schoenberner Sz Napiwotzki 
1990). Another six were rejected for various reasons related to contamination of their photometry: 

• WD 0457-103 is in a known spectroscopic binary with a bright G-K star (63 Eri; Shara et al. 
1997; Barstow et al. 1994); it remains unresolved despite a high resolution Hubble Space 
Telescope imaging investigation (Barstow et al. 2001). 
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• WD 0725+318 has a nearby source (RA=07:28:11.65, dec=+31:43:46.57, J2000), which is 
revealed in the SDSS images of the field as a red background galaxy that is very bright in the 
IR and likely contaminates the WD photometry in both 2MASS and WISE. 

• WD 1109-225 is in a known unresolved spectroscopic binary with a bright A star (f3 Crt; 
Fleming et al. 1991). 

• WD 1201+437 is classified as a DC+dMe binary by Fleming et al. (1993). Xu et al. (1999) 
classify it as a quasar based on its X-ray properties; however, we note that their X-ray error 
circle has a 9" radius, so this could be a mis-identification. In either case, we would remove 
it from our sample. 

• WD 1235+321 has a stellar profile in the WISE images that is faint and extended, while 
the SDSS images show two nearby faint point sources that are likely contaminating the 
WISE photometry. These sources are: a star located 5.6" west with i = 20.11 mag and 
z = 19.94 mag, and a galaxy located 6.2" east with i = 19.87 mag and z = 19.64 mag. There 
are three additional faint sources within 13" of the WD (two to the south, one to the southwest 
in the direction of the bright star located 25" southwest of the WD), all with i > 21 mag and 
z > 22 mag. 

• WD 1859+429 is in a wide binary with a common proper motion companion located « 15" 
to the northeast. The common proper motion companion is not problematic; however, due 
to its proper motion, the WD is superposed on a field star in the 2MASS images. It is likely 
that this star contaminates the WISE photometry. 

This leaves 10 WDs that we classified as targets-of-interest and for which we constructed UV- 
IR spectral energy distributions (SEDs). In addition to the 2MASS (Skrutskie et al. 2006) All Sky 
Data Release Point Source Catalog and WISE All Sky Release 2 photometry, we utilized photometric 
data from the Galaxy Evolution Explorer ( GALEX ; Martin et al. 2005), the Sloan Digital Sky 
Survey (SDSS Data Release 7; Abazajian et al. 2009), the AAVSO Photometric All Sky Survey 
(APASS Data Release 6) 3 , and the Spitzer Space Telescope Enhanced Imaging Products Source List 
(Cryogenic Release v2.0, January 2013) 4 . For the purpose of modeling the SEDs, we converted the 
various photometric measurements from magnitudes into flux densities using published zero points 
for each photometric band (Bessell 1979; Skrutskie et al. 2006; Abazajian et al. 2009; Wright et al. 
2010). All of the photometric data for the targets-of-interest plus the CSPNs and known WDs with 
dust disks, along with published spectral types, effective temperatures, and surface gravities of the 
WDs, are listed in Tables 1-3. 


2 See the Explanatory Supplement at http://wise2.ipac.caltech.edu/docs/release/allsky/expsup/. 

3 See http://www.aavso.org/apass. 

4 See the Explanatory Supplement at 

http : //irsa. ipac . caltech. edu/data/SPITZER/Enhanced/Imaging/docs/Spitzer_EIP_expsup .pdf . 
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We then fit either a DA or DB model, as appropriate for the published WD type, to the UV- 
optical-near-IR ( JHK) portion of the target SEDs, using grids of H- and He-rich WD cooling models 
(kindly provided by P. Bergeron) that include the GALEX, SDSS, Johnson UBVRI, 2MASS, and 
WISE bands (Bergeron et al. 1995; Holberg &; Bergeron 2006). The published WD temperature 
and log g values from Table 1 were used as initial values (a “typical” value of log g = 8.0 was assumed 
in cases for which there is no published value). We searched for best- fitting models within the lcr 
uncertainties of the literature values of T e g- and log g. For WDs with no published uncertainty for 
T e ff, we assumed ±2000 K. 

Seven of the targets have an obvious excess over the WD model in the IR. For these objects, 
we re-fit the entire SED using an additional circumstellar dust disk component. The model dust 
disk SED was calculated as originally devised by Jura (2003), following the procedure described 
in Debes et al. (2011), with free parameters of inner radius, width, and inclination. The minimum 
allowed inner disk radius was given by a conservative sublimation temperature of T su bi = 2100 K; 
in the case of our hottest WD (WD 0420±520, T e g = 24300 K), we relaxed this criterion to 
Tsubi = 2500 K in order to obtain the best dust disk model fit 5 . 

As noted above, the broad 10 /im silicate emission feature that is commonly seen in the 
IR spectra of WDs with circumstellar dust (Jura et al. 2009) could contribute in the W3 band; 
consequently, we did not use the W3 photometry to constrain the models. Any model that is too 
faint at W3 indicates the presence of a silicate emission feature or a significant amount of cool 
dust at large distances from the WD. If the model is too bright at W3, then this likely indicates 
that the assumed outer radius of the disk is too large. The silicate emission feature at 18-20 n nr 
(e.g., seen in the IR spectrum of the archetype dusty WD G29-38 = WD 2326±049; Reach et al. 
2009) is generally weaker in amplitude but broader than the 10 /jrn feature (Papoular & Pegourie 
1983; Thompson et al. 2003; Jura et al. 2007b) 6 , and could contribute in the W4 band. The disk 
inclinations are generally poorly constrained by the models. The SEDs and model fits are shown 
in Figure 2; the WD and dust disk model parameters are listed in Table 4. These seven objects are 
our new candidate WDs with dust disks. 

There are three reduced chi-squared (x 2 ) values listed in Table 4 for each model. These provide 
different measurements of the goodness of the model fit: 

• The parameter x(/ d re f ers to the goodness of just the WD model component compared to 
only the UV-optical-near-IR (JHK) data. In this wavelength region, we expect generally 


5 Silicate dust, for example, is generally assumed to sublimate at T > 1500-2000 K (e.g., Pollack et al. 1994; 
Kobayashi et al. 2011); however, Rafikov & Garmilla (2012) recently calculated that due to the high metal vapor 
pressure at the inner edge of WD dust disks, the dust grain sublimation temperatures can be 300-400 K higher than 
is generally assumed, and the inner edge of the disk can be superheated to as high as 2500-3500 K. 

6 Laboratory experiments suggest that “fresh” silicates have a very large 10/20 /im flux ratio, but factors such as 
metal content or oxidation - that is, age - of the silicate compounds can increase the strength of the 20 /im feature 
relative to the 10 /im feature (Nuth & Hecht 1990). Dust temperature can also affect this ratio (Suh 1999). 
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good agreement with a “naked” WD model even in the presence of circumstellar dust (which 
contributes most strongly at mid-IR and longer wavelengths). Large values of this statistic 
likely indicate deviations from the WD model in the UV; the GALEX data are typically 
among the brightest points in the SED and have very small relative errors. If x^ d l ar § e 
because of deviations in the UV, then the values of the other two x 2 statistics will also be 
large, as the errors are dominated by the poor fit in the UV. For example, removing the two 
UV points from the WD 2329+407 SED reduces x^ d from 107 to 2.1 and x disk ( see below) 
from 94 to 2.6. In the case of WD 1046-017, removing the UV points causes all of its x 2 
values to drop to ~ 1.4. 

• The parameter x 2 u refers to the goodness of just the WD model component compared to all 
of the available photometric data. Large values (e.g., x 2 u > Xwd > 1) indicate the need for 
an additional model component in the IR. 

• The parameter Xdisk re f ers to the goodness of the WD + dust disk model compared to all 
of the available photometric data. If Xdisk < Xain then the model was improved by the 
addition of a dust disk component. As noted above, this statistic can be large because it is 
dominated by a poor fit in the UV. It can also be large because the W3 (and/or W4) data 
points are bright compared to the model (see above, as well as individual target notes below); 
for example, removing the bright W3 and W4 points from the SED of WD 0420+520 reduces 

Xdisk from 7 - 8 to °- 8 - 


Three targets-of-interest selected by our criteria (WD 1146-290, WD 1330+473, and WD 2152- 
548) show no strong IR excess (see Figure 3), and they are among the selected WDs with (Wl— W2) 
color closest to the +0.3 mag criterion (see Figure 1). The first of these is a very cool WD (see 
Table 1 and notes below). It has a distinctive SED shape compared to the others in our selected 
sample, and its selection by our criteria was likely a by-product of its very low temperature. The 
latter two WDs have only very slightly elevated W2 flux densities compared to the WD model. In 
particular, only the W2 point in WD 2152-548 is elevated compared to the model; the adjacent 
IRAC-2 point (as well as the other IRAC data) agree with the model. Thus, the W2 value for 
WD 2152-548 should be treated with caution. If real, then the IR excesses in WD 1330+473 and 
WD 2152-548 are very weak and could indicate that only a small amount of cool dust is present. 
This situation presents similarities to the weak IR excesses found in PG 1457-086 (Farihi et al. 
2009) and HE 0106-3253 (Farihi et al. 2010), which are inferred to be due to narrow circumstellar 
rings of dust instead of full disks, as well as to the several known WDs with combined gas and dust 
disks (e.g., Gansicke et al. 2006; Brinkworth et al. 2009, 2012). The presence of co-mingled gas and 
dust disks could indicate that a significant amount of dust has been converted to gas through either 
sublimation (Melis et al. 2010) or collisions (with itself or possibly with pre-existing circumstellar 
material; Jura 2008; Girven et al. 2012). 


2.1. Notes on Individual Targets 


We have examined the publication record for each of the ten targets-of-interest, and briefly 
describe any relevant features below. We also vetted the WISE data for each of these sources for 
evidence of red contaminating sources in the photometry aperture (as described in Debes et al. 
2011); relevant notes are included below. 


2.1.1. WD 0249-052 

Voss et al. (2007) and Limoges & Bergeron (2010) do not note any atmospheric contamination 
(including hydrogen) in this DB WD. In the former study, hydrogen was assumed to be absent 
unless a visual inspection of the optical spectrum revealed H lines (corresponding to a detection 
limit of Ha equivalent width > 300 mA). In the latter study, it appears that hydrogen was only 
utilized in the model spectrum analysis if the WD had been previously identified as a hydrogen- 
rich helium (DBA) WD. Additional examination of the two (somewhat noisy) spectra of this WD 
from the Voss et al. (2007) study yields no obvious metal lines, and limits of [Ca/He] < —8.0 and 
[Mg/He] < —6.7 (D. Koester, private communication). On the other hand, Bergeron et al. (2011) 
find [H/He]= —5.47(59) (but no metals) by utilizing high S/N spectra of the Ha region. None of 
these three analyses, however, conclude that metal contamination (signified by the presence of Ca 
absorption in the optical spectrum) is present. We have no concerns about the quality of the WISE 
photometry. 


2.1.2. WD 0420-731 

There is no detailed information on this WD in the literature. However, there is a source 
(WISE-J041933. 70-730333. 9) located ~ 22" northwest of the WD, which is faint in W1 and W2, 
but becomes much brighter in the W3 and W4 bands (slightly brighter than the WD). This source 
separation is well beyond the 1.3 x (FWHMwi) ~ 7.8" radius, interior to which contamination of the 
WISE All Sky Catalog photometry can occur (see discussion in Debes et al. 2011); nonetheless, we 
tested for contamination from the neighboring source. To do so, we used the IRAF 7 implementation 
of DAOPHOT (Stetson 1987) to obtain PSF-ht photometry for the WD and the nearby source in the 
W3 Atlas images, using the nearby bright star WISE-J041948. 50-730317.2 as a PSF template and 
magnitude calibrator. We obtain W3 ps f = 11.69(40) mag for the WD and W3 ps f = 11.60(33) mag 
for the nearby source, in agreement with the WISE catalog values of W3 = 11.700(112) mag and 
W3 = 11.642(106) mag, respectively. Nonetheless, the local background is patchy and bright in 
W3 and W4, so it is prudent to treat the W3 and W4 photometry as upper limits until higher 


'IR.AF is distributed by the National Optical Astronomy Observatories, which are operated by the Association of 
Universities for Research in Astronomy, Inc., under cooperative agreement with the National Science Foundation. 
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resolution imaging data are available. 


2.1.3. WD 0420+520 

There is no detailed information on this WD in the literature. The cataloged W3 and W4 flux 
densities for this target are quite bright, and there is no obvious point source at the position of the 
WD in the WISE W3 and W4 Atlas images. Thus, these values should be treated as upper limits. 


2.I.4. WD 0836+404 

This is a ZZ Ceti type pulsating WD (Vauclair et al. 1997). Farihi et al. (2005) found no 
evidence for a low luminosity companion from a survey utilizing proper motion measurements, 
deep imaging, and near-IR photometry. A limit on atmospheric metal contamination was set by 
Zuckerman et al. (2003), at [Ca/H]< —7.72. 

There is a bright (V=10.9 mag) star (2MASS J08401164+4015211) located « 43" east of 
the WD. While this star is far enough from the WD to not pose a contamination risk for the 
WISE photometry, we note that a diffraction spike from the star passes near the WD in the 
WISE Atlas images. Contamination warnings due to diffraction spikes are included in the WISE 
All Sky Catalog; such a warning was not flagged for WD 0836+404. Nonetheless, to confirm 
this we performed DAOPHOT PSF-subtraction photometry on the W2 Atlas image, as described 
above. We used the nearby stars WISE-J084000. 92+401704.4 and WISE-J084022. 86+401250. 6 
to construct a PSF template, and the mean photometry of the PSF stars plus several nearby 
stars ( WISE- J083954. 52+401509. 5, WISE-J084001.71+401415.7, and WISE-J084022.35+401424.8) 
that are comparable in brightness to the WD as a magnitude calibrator. We obtain W2 ps f = 
15.29(51) mag for the WD, in agreement with its WISE catalog value of W2= 15.245(115) mag. 
The nominal W2 PSF photometry for all 6 of the tested stars (including the WD) agrees to better 
than 1% with the WISE catalog values. So, we consider it unlikely that the nearby diffraction spike 
has contaminated the WISE photometry of the WD. 


2.1.5. WD 1046-017 

This is a known DB WD and there is tentative evidence that it might be nretal-rich: Sion et al. 
(1988) noted a possible weak Ca II K feature in its optical spectrum (equivalent width < 10 mA), 
while Zuckerman et al. (2010), Bergeron et al. (2011), and Jura & Xu (2012) set limits of [Ca/He] 
< —10.9 (Ca II equivalent width < 9 mA), [H/He] < —6.5, and log(dM meta \ s / dt [g s^ 1 ]) < 6.20. 
Thus, the metal-rich status of this WD remains uncertain (but unlikely). Farihi et al. (2005) found 
no evidence for a low luminosity companion from a survey utilizing proper motion measurements, 


deep imaging, and near-IR photometry. We have no concerns about the quality of the WISE 
photometry. 


2.1.6. WD 1146-290 

The equivalent width of Ha in this DA WD is 5.9 A (Bergeron et al. 1997). There is no other 
detailed information on this WD in the literature. We have no concerns about the quality of the 
WISE photometry. 


2.1.7. WD 1330+473 

No IR excess or evidence for a dust disk is noted in the near- and mid-IR photometric and spec- 
troscopic survey of Barber et al. (2012). Farihi et al. (2005) found no evidence for a low luminosity 
companion. We have no concerns about the quality of the WISE photometry. 


2.1.8. WD 1448+411 

There is no detailed information on this WD in the literature. We have no concerns about the 
quality of the WISE photometry. 


2.1.9. WD 2152-548 

This object was first discovered as an X-ray source in the Einstein satellite slew survey 
(Elvis et al. 1992), and was later confirmed as a hot DA WD optical counterpart to a ROSAT 
X-ray source (Mason et al. 1995; Marsh et al. 1997). Bannister et al. (2003) report possible weak 
photospheric metal contamination in this WD from high resolution UV spectroscopic observations, 
and note that it is “an object deserving of further attention.” Dickinson et al. (2012) re-examined 
the origin of previously reported circumstellar features in 23 hot DA WDs. While unambiguous 
re-detections of circumstellar material were made for eight other WDs, they were unable to confirm 
the Ca II contamination reported by Bannister et al. (2003) for WD 2152-548. We have no concerns 
about the quality of the WISE photometry. 


2.1.10. WD 2329+407 

This WD was noted as non-magnetic in the spectropolarimetric survey of Schmidt Sz Smith 
(1995). Farihi et al. (2005) found no evidence for a low luminosity companion. There is a faint field 
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star located ~ 8" southwest of the WD. Proper motion of the WD between the DSS-1 and -2 epochs 
shows that this object is not related to the WD. The neighboring source is visible in the 2MASS im- 
ages but is not listed in the 2MASS Point Source Catalog. The W1 Atlas image shows a slight exten- 
sion to the WD PSF that is consistent with the presence of this source, but is not visible in the W2 
and W3 Atlas images (the WD is not visible at all in the W4 Atlas image; the W4 photometry for the 
WD is an upper limit). As with WD 0420-731, we performed DAOPHOT PSF-subtraction photom- 
etry, this time on the W1 Atlas image. We used the nearby bright star WISE- J233143. 17+410133.1 
as a PSF template, and the mean photometry of several nearby stars (WISE-J233140. 08+410100. 9, 
WISE- J233134. 97+410217.7, WISE-J233129.20+410124.4, and WISE-J233135.68+410012.5) that 
are comparable in brightness to the WD as a magnitude calibrator. The corresponding PSF- 
subtraction of the WD was successful (see Figure 4), and we obtain Wl ps f = 13.75(5) mag, which 
is in agreement with the catalog photometry of Wl= 13.757(27) mag. Thus, there appears to be 
no overt reason to be concerned about the WISE photometry of WD 2329+407. 


3. Discussion and Conclusions 

Many of the “original” WDs with dust disks that were discovered before 2011 have IR excesses 
that are detectable in the near-IR (JHK) bands. In all but one (WD 0420+520) of the new dust disk 
candidates presented here, the excess emission due to dust is only apparent at wavelengths > 3 /an. 
It is possible that the current census of WDs with dust disks that produce an excess detectable 
at K-band and shorter wavelengths (e.g., using 2MASS or UKIDSS data; see Wachter et al. 2003; 
Wellhouse et al. 2005; Hoard et al. 2007; Steele et al. 2011) is close to complete for the entire sample 
of known WDs (at least to the detection limits of existing near-IR all-sky surveys) . The WD dust 
disk candidates now being found using longer wavelength data from WISE and Spitzer are drawn 
from a previously underrepresented region of parameter space, in which the dust disks are overall 
cooler, narrower in radial extent, and/or contain fewer emitting grains. 

In Figure 5, we have plotted the dust disk inner edge radius as a function of WD effective 
temperature for the seven new candidates from this paper, the candidates in Debes et al. (2011), 
and various published dust disk models for other WDs (see Table 5). There is a direct relationship 
between the WD temperature and the temperature of dust at a given radius (Jura 2003), 

T dust (R) ocir 3 / 4 T wd , (1) 

which is used to plot isotemperature contours in Figure 5 for the model dust disks. For an assumed 
dust sublimation temperature, the corresponding contour shows the minimum inner radius of the 
disk as a function of WD temperature; in general, the contours illustrate the temperature at a 
given radius in the disk for a given WD temperature. A typical boundary condition of WD dust 
disk models is that no dust is allowed to be hotter than the assumed sublimation temperature; in 
other words, dust is not allowed at radii closer to the WD than the radius at which the sublimation 
temperature (according to equation 1) is reached. In many cases, the hottest dust in the disk 
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models is at the assumed sublimation temperature (disk models with inner edges at the assumed 
sublimation temperature are indicated in Tables 4 and 5). Such a disk extends inward toward the 
WD as far as possible - its inner edge lies at the sublimation radius around the WD, so dust cannot 
survive any closer to the WD. On the other hand, disks that have inner edge radii corresponding 
to temperatures below the dust sublimation temperature could, in principle, extend inward closer 
to the WD but do not. In such cases, possible reasons for the lack of dust grains close to the WD 
include (but are not limited to): 

• Grains close to the WD might have been depleted due to a higher rate of collisions with other 
grains. This might have occurred in the handful of known WDs with gas+dust disks, in which 
the gas and dust share spatially overlapping, but not identical, radial distributions, implying 
that dust sublimation alone cannot account for the presence of gas/lack of close-in dust (e.g., 
Brinkworth et al. 2012). Additionally, objects like WD 1456+298 (G166-158) have a very 
weak IR excess corresponding to the presence of only cool dust (see additional discussion of 
this object below). 

• A spinning WD with a strong magnetic field might sweep up paramagnetic or diamagnetic 
dust grains interior to a critical radius. A similar process operates in the intermediate polar 
class of cataclysmic variable to produce a truncated gaseous accretion disk around the WD 
(see review of this class in Warner 2003). Based on observations of the Ca II lines during 
pulsations of the archetype dusty WD G29-38, Thompson et al. (2010) have suggested that 
the Ca is being accreted onto the poles of the WD, rather than equatorially, suggesting that 
the WD is magnetic. Brinkworth et al. (2007) (see the Appendix in that paper) calculated 
the critical surface charge on a dust grain, Q C rit, such that the motion of the dust would be 
influenced by the WD magnetic field. They found that Q C rit for dust grains at the sublimation 
radius around a 14,200 K WD with a field strength of B ~ 25 MG is more than an order 
of magnitude larger than the likely surface charge of the grains based on observations of in 
situ interplanetary dust grains in the Solar System. However, they note that the value of 
Q cr it is several orders of magnitude smaller than the value £") max , at which the electrostatic 
tensile stress in the dust grain interiors would be sufficient to fracture (i.e. , destroy) the 
grains. In addition, Q C rit oc R 3 / 2 B (where R. is the distance from the center of the WD 
to the dust grain), leaving open the possibility that magnetic interactions could be effective 
around WDs with larger magnetic fields and/or cooler temperatures (allowing the dust to 
approach closer to the WD before sublimating). Finally, Brinkworth et al. (2007) did not 
consider any intrinsic para- or diamagnetic properties of the dust itself, which could enhance 
interaction with the WD magnetic field. Other than these two examples (one observational, 
one theoretical), there has been (to our knowledge) little exploration of this possibility. 

There is, of course, some ambiguity in the interpretation of Figure 5. A given disk could appear 
to not reach the sublimation radius if the assumed sublimation temperature was higher than the 
true sublimation temperature for the particular species of dust in the observed disk. In the bulk of 
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cases, however, the published models of dust disks extend inward to the sublimation radius even 
when the assumed sublimation temperature is quite high (e.g., a majority of the large sample in 
Debes et al. 2011 with T su ^i = 2100 K, or WD 0420+520 in this work). Apparently, most of the 
currently known dusty WDs have “hot” disks in which the dust extends inward quite close to the 
WD, until it reaches the “sublimation barrier.” In contrast, five of the seven dust disk candidates 
presented here have maximum dust temperatures of < 1000 K, safely below reasonable lower limits 
of the sublimation temperature for metallic dust. They are, therefore, truly “cool” disks in which 
dust is depleted close to the WD and is mainly present substantially exterior to the sublimation 
radius. These disks are narrower in radial extent than an otherwise identical disk in which the inner 
edge extends all the way to the sublimation radius. The apparent bias toward sublimation- limited 
disks does not, however, imply that “hot” WD dust disks are necessarily more intrinsically common 
than “cool” disks. The “hot” disks are easier to find, since they produce strong IR excesses that 
can be detected in the near-IR. 

Three of our dust disk candidates (WD 0249-052, WD 0836+404, WD 1046-017) have published 
optical spectroscopic studies noting the absence of atmospheric metal pollution. One would expect 
to find metals in the atmosphere of a WD with an IR excess indicating a dust disk, as some of the 
dust will accrete onto the star. The DB WDs (WD 0249-052, WD 1046-017) have gravitational 
settling times for metals in their helium-rich atmospheres that are much longer than for DA WDs; 
metals can persist for up to a Myr or longer after accretion ceases, making them appear as metal-rich 
WDs for a substantial time after the dust disk has dissipated (typical disk lifetimes are ~0.03-5 Myr; 
Girven et al. 2012). While there appear to be firm constraints on the lack of accreted metals in 
WD 0249-052 (see Section 2.1.1), the situation for WD 1046-017 is less certain (see Section 2.1.5). 
Possibly the most viable explanation for the presence of a dust disk around a DB WD that is not 
metal-polluted is that the disk is newly formed and the abundance of accreted metals on the WD 
is not yet in a steady state and has not reached the threshold of detection. 

On the other hand, because of the short settling times of DA WDs, essentially any process 
that causes the accretion of metals from the dust disk to be out of steady state could produce 
the phenomenon of a WD with an IR excess but no spectroscopic signature of metal pollution. 
In the case of the DA WD among these three, we note that the limit on the metal contamina- 
tion of WD 0836+404 is not very stringent ([Ca/H]< —7.7; Zuckerman et al. 2003). This leaves 
some ambiguity in the conclusion that the WD has no metal pollution; for example, the SED 
of WD 0836+404 (like the other two; also see discussion above) is similar in appearance to that 
of G166-158. The latter is a metal-polluted DA WD with a very weak IR excess that only be- 
comes apparent longward of ~ 5 //m, and likely indicates a narrow annulus of cool dust located 
relatively far from the WD with a large, inner, dust-free zone. Zuckerman et al. (2003) measured 
[Ca/H]= —9.3 for G166-158, implying that WD 0836+404 could still be considered metal-rich with 
[Ca/H]< —7.7, even if the true Ca abundance limit is an order of magnitude lower, given that 
Zuckerman et al. (2003) assumed a somewhat higher temperature for WD 0836+404. By the same 
token, WD 0836+404 (~ 11700 K) is hotter than G166-158 (~ 7400 K), so the former would require 
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a higher Ca abundance to show lines than the latter (i.e., they could have comparable levels of Ca 
enrichment, but only the cooler WD might show Ca lines in its optical spectrum). Regardless, a 
DA WD with a dust disk but no metal pollution might be explained by a number of scenarios, 
including the inner disk dust depletion and “magnetic sweeping” scenarios discussed above, as well 
as non-steady state accretion from a newly-formed dust disk (as suggested above for the case of 
DB WDs). 

We end with a final note of caution. As discussed in Debes et al. (2011), the large WISE 
PSF can defy even the most careful vetting process and allow unidentified contamination from 
faint, red, unresolved sources. Thus, the seven WDs highlighted here should be considered only as 
dust disk candidates until independent confirmation of their IR excesses can be obtained. This is 
especially true for the three WDs, discussed above, that do not appear to show the signature of 
metal contamination in their optical spectra. In these cases, the existing IR data do not strongly 
constrain a dust disk model. Additional longer wavelength mid-IR observations are needed to rule 
out the possibility that a very cool brown dwarf companion could be responsible for the observed 
IR excess. 

This work is based on data, data products, and other resources obtained from: (a) The Two Mi- 
cron All Sky Survey (2MASS), a joint project of the University of Massachusetts and the Infrared 
Processing and Analysis Center (IPAC) /California Institute of Technology (Caltech), funded by 
the National Aeronautics and Space Administration (NASA) and the National Science Foundation 
(NSF). (b) The American Association of Variable Star Observers (AAVSO) Photometric All-Sky 
Survey (APASS), funded by the Robert Martin Ayers Sciences Fund, (c) NASA’s Astrophysics 
Data System, (d) The Galaxy Evolution Explorer ( GALEX ), a NASA Small Explorer launched in 
April 2003 and operated for NASA by Caltech under NASA contract NAS-98034. The GALEX 
data products were obtained from the Multimission Archive at the Space Telescope Science In- 
stitute (MAST). STScI is operated by the Association of Universities for Research in Astronomy, 
Inc., under NASA contract NAS5-26555. Support for MAST for non-HST data is provided by the 
NASA Office of Space Science via grant NNX09AF08G and by other grants and contracts, (e) The 
NASA/IPAC Infrared Science Archive (IRS A), which is operated by the Jet Propulsion Laboratory 
(JPL), Caltech, under a contract with NASA, including the Spitzer Enhanced Imaging Products 
based on observations obtained with the Spitzer Space Telescope, which is operated by JPL, Cal- 
tech, under a contract with NASA, (f) The SIMBAD database, operated at CDS, Strasbourg, 
France, (g) The Sloan Digital Sky Survey (SDSS and SDSS-II), whose funding has been provided 
by the Alfred P. Sloan Foundation, the Participating Institutions, the NSF, the U.S. Department of 
Energy, NASA, the Japanese Monbukagakusho, the Max Planck Society, and the Higher Education 
Funding Council for England. The SDSS Web site is http:/ /www. sdss.org/. The SDSS is managed 
by the Astrophysical Research Consortium for the Participating Institutions. The Participating 
Institutions are the American Museum of Natural History, Astrophysical Institute Potsdam, Uni- 
versity of Basel, University of Cambridge, Case Western Reserve University, University of Chicago, 
Drexel University, Fermilab, the Institute for Advanced Study, the Japan Participation Group, 
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Johns Hopkins University, the Joint Institute for Nuclear Astrophysics, the Kavli Institute for Par- 
ticle Astrophysics and Cosmology, the Korean Scientist Group, the Chinese Academy of Sciences 
(LAMOST), Los Alamos National Laboratory, the Max-Planck-Institute for Astronomy (MPIA), 
the Max-Planck- Institute for Astrophysics (MPA), New Mexico State University, Ohio State Uni- 
versity, University of Pittsburgh, University of Portsmouth, Princeton University, the United States 
Naval Observatory, and the University of Washington, (h) The Wide-field, Infrared Survey Explorer 
(WISE), which is a joint project of the University of California, Los Angeles, and JPL, Caltech, 
funded by NASA. 
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Fig. 1. Infrared color-color diagram of MS99 WDs with 2MASS J and WISE W1 and W2 
detections. The vertical dashed line marks a (Wl— W2) color of +0.3 mag, which was used in the 
target selection process. Representative color index error bars are shown in the upper left; from 
top to bottom, 25%, 50%, 75%, and 90% of the detected targets have photometric uncertainties 
smaller than the indicated error bars. The points are symbol-coded as follows: new candidate WDs 
with dust disks (blue filled circles), selected targets that are not dust disk candidates (blue unfilled 
circles), known WDs with dust disks that are in (green squares) or not in (black squares) MS99, 
known unresolved binaries (red diamonds), central stars of planetary nebulae (red downward facing 
triangles), WISE photometry is contaminated (red upward facing triangles), and remaining WDs 
that did not satisfy our selection criteria as targets-of-interest (small gray circles). WD 1201+437 
is also a known unresolved binary with (J— Wl) = +3.8 mag and, for clarity, is the only target not 
plotted here. 
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Fig. 2. — Spectral energy distributions of the seven new WD with dust disk candidates. Photometric 
values are shown as filled circles; cataloged 2a (95% confidence) upper limits for non-detections are 
shown as downward arrows. The unfilled circles are W3 and W4 photometry that is of questionable 
quality and should be treated as upper limits (see Sections 2.1.2 and 2.1.3). The wavelengths of the 
various photometric bands are indicated at the top of the figure (see Tables 1-3 for the wavelength 
of each band). The dashed line is a WD model fit, while the solid line is a combined WD + dust disk 
model (see Table 4 for disk parameters). The two models are generally indistinguishable shortward 
of « 2 //m. 
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Fig. 3. — As in Figure 2, but for the three WDs selected by our target criteria that are not dust 
disk candidates. 



Fig. 4. — The WISE W1 Atlas image of WD 2329+407 before (left panel) and after (right panel) 
PSF-subtraction of the WD. The circle is centered on the WD and shows the extent of the PSF 
radial profile fit. The image is ~ SO" in both dimensions, with a plate scale of « 1.4" pixel -1 , and 
is oriented north up, west to the right. This figure is representative of all of the PSF subtraction 
photometry performed here, as discussed in Section 2.1. 
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Fig. 5. — Inner edge radius as a function of WD effective temperature for published WD dust disk 
models. The dashed lines are isotemperature contours that show the temperature at a given radius 
in the disk for a given WD temperature. Models that are “sublimation- limited” (i.e. , the dust 
located closest to the WD is constrained to be no hotter than the sublimation temperature) will he 
along the isotemperature contour corresponding to the dust sublimation temperature assumed in 
the model. The small filled circles are the new WD dust disk candidates from Debes et al. (2011) 
(with assumed sublimation temperature of 2100 K, shown as a dotted isotemperature contour) and 
the large filled circles are the new candidates from this work. The other points are previously known 
WD dust disks, using published disk models and WD effective temperatures (see Table 5; some 
WDs are plotted more than once corresponding to different published dust disk model parameters) , 
and are symbol-coded as follows: unfilled triangles represent models with assumed sublimation 
temperatures of 1200-1500 K, unfilled squares have assumed sublimation temperature of 1800 K 
(shown as a dotted isotemperature contour), and unfilled circles have various assumed sublimation 
temperatures in the range 1000-2000 K. Disk radii were rounded to integers, so some plotted 
points fall at radius values slightly smaller or larger than the minimum radius for the corresponding 
assumed sublimation temperature. See text for additional discussion. 
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Table 4. Model Parameters for New WD Dust Disk Candidates 


WD 

Other Name 

^eff,wd 

(K) 

distance 

(pc) 

-^in,disk 

(-ftwd) 

-^out,disk 

(-ftwd) 

inclination 

n 

Xwd a 

V 2 b 

-^all 

xLc 

0249-052 

HE 0245-0514 

17823 

104 

30 

52 

80 

1.3 

5.1 

1.9 

0420-731 


17653 

79 

18 

88 

71 

2.8 

151 

6.1 d 

0420+520 

KPD 0420+5203 

24301 

76 

12* 

28 

71 

1.4 

64 

7.8 d 

0836+404 

DF Lyn 

11712 

59 

16 

86 

89 

1.0 

2.7 

1.2 

1046-017 

GD 124 

14266 

75 

28 

38 

80 

17. 6 e 

18. 0 e 

17. 7 e 

1448+411 

CBS 204 

13571 

80 

26 

36 

80 

0.3 

9.5 

0.9 d 

2329+407 

EGGR 160 

13900 

33 

28 

68 

80 

107 e 

155 e 

94 e 

1146-290 

Ruiz 440-146 

5000 

26 




12 

32 


1330+473 

PG 

21223 

91 




461 

384 


2152-548 

1ES 2152-54.8 

45050 

123 




0.7 

10 



* Dust located at R[ n disk is at the assumed sublimation temperature. 

a x 2 value of the WD model fit compared to only the GALEX UV, optical, and 2MASS near-IR data points. 
b x 2 value of the WD model fit compared to all available data points. 
c x 2 value of the WD + dust disk model fit compared to all available data points. 

d These values are large due to the inclusion of bright W3 and/or W4 points — see discussion in the text. 
e Large x 2 values due to poor fit in the UV - see discussion in the text. 
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